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Abstract 

It has been found that recent results on forward jet production from deep 
inelastic scattering can neither be reproduced by models which are based on 
leading order a s QCD matrix elements and parton showers nor by next-to- 
leading order calculations. The measurement of forward jet cross sections has 
been suggested as a promising probe of new small x parton dynamics and the 
question is whether these data provide an indication of this. The same question 
arises for other experimental data in deep inelastic scattering at small x which 
can not be described by conventional models for deep inelastic scattering. In 
this paper the influence of resolved photon processes has been investigated and 
it has been studied to what extent the inclusion of such processes in addition 
to normal deep inelastic scattering leads to agreement with data. It is shown 
that two DGLAP evolution chains from the hard scattering process towards the 
proton and the photon, respectively, are sufficient to describe effects, observed 
in the HERA data, which have been attributed to BFKL dynamics. 



1 Introduction 

Experimental data from deep inelastic scattering (DIS) in a kinematic region where new 
parton dynamics is expected to become noticeable i.e. at small values of the scaled proton 
momentum, x, are not described by models based on interactions with pointlike photons. 
In a previous paper Jl| we have demonstrated that the addition of interactions through 
resolved photons offers a possible explanation to the observed discrepancies and leads to 
good agreement with all available data. 

This paper is devoted to a more detailed discussion of the resolved photon concept 
and comparisons with data on forward jet production in DIS since the forward jet cross 
section has been advocated as a particularly sensitive measure of small x parton dynamics 
|2|, [|. Analytic calculations based on the BFKL equation in the leading logarithmic 
approximation (LLA) are in fair agreement with data. However, recent calculations of 
the BFKL kernel in the next-to-leading logarithmic approximation (NLLA) M have given 
surprisingly large corrections, and it remains to be shown whether the data can still be 
reasonably described. 

Monte Carlo generators based on direct, point-like photon interactions (DIR model), 
calculated from leading order (order a s ) QCD matrix elements, and leading log parton 



showers based on the DGLAP evolution do not take any new parton dynamics in the 
small x region into account and are therefore not expected to fit the experimental data. 
Recent results from the HI || and ZEUS pi experiments on forward jet production 
exhibit significant deviations from the predictions of such models. Also next-to-leading 
order calculations (NLO i.e. order a 2 ) predict too small a cross section compared to data. 

The study of forward jet production with contributions from direct as well as resolved 
photon processes has been performed using the RAPGAP 2.06 J7], ^ Monte Carlo event 
generator. 



2 Resolved Photons in DIS 

In electron-proton scattering the internal structure of the proton as well as of the ex- 
changed photon can be resolved provided the scale of the hard subprocess is larger than 
the inverse radius of the proton, l/R 2 ~ A.q CD , and the photon, l/R 2 ~ Q 2 , respectively. 
Resolved photon processes play an important role in photo-production of high jets, 
where Q 2 ~ 0, but they can also give considerable contributions to DIS processes |9|, [L(| 
if the scale /i 2 of the hard subprocess is larger than Q 2 , the inverse size of the photon. 

In the following we give a brief description of the model for resolved virtual photons 
used in the Monte Carlo generator RAPGAP. Given the fractional momentum transfer 
of the incoming electron to the exchanged photon, the Equivalent Photon Approximation 
provides the flux of virtual transversely polarized photons || 0, and references therein]. 
The contribution from longitudinally polarized photons has been neglected. The partonic 
structure of the virtual photon is defined by parameterizations of the parton densities, 
x^/ 7 (x 7 , /i 2 , Q 2 ), which depend on the two scales /i 2 and Q 2 ||TT], [1^, The following 
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hard subprocesses are considered (RES model): gg — > qq, gg — > gg, qg — > qg, qq — > gg, 
qq — > qq, qq — > qq. Parton showers on both the proton and the photon side are included. 
The generic diagram for the process q^g v — * qg including parton showers is shown in 
Fig. I- 

Since the photon structure function depends on the scale, fi 2 , of the hard scattering 
process, the cross section of resolved photon processes will consequently also depend on 
the choice of this scale. It has to be carefully considered in which range of \j?IQ\ the 
photon-parton cross section can be factorised into a parton-parton cross section convoluted 
with the parton density of the photon. The parton density of the photon is evolved 
from a starting scale Q\ to the scale /x 2 , the virtuality at the hard subprocess, giving a 
resummation to all orders. 



2.1 Parton Distribution Functions 

Due to factorization of the cross sections the parton densities of both the virtual photon 
and the proton enter into the calculations. The proton structure function, F-i, has been 
measured to high accuracy and therefore the various parameterizations only give marginal 
differences in the measurable kinematic region. Two parameterizations of the parton 
distribution in the proton, GRV 94 HO (DIS) and CTEQ4D have been considered, which 
both give good agreement with the proton structure function data [13, IS|. It was found 




Figure 1: Deep inelastic scattering with a resolved virtual photon and the q^g p — > qg 
partonic subprocess. 



that the produced results were identical within the percent level, when keeping Aqcd 
fixed. In the following we use only GRV HO. 

The photon can interact via its partons either in a bound vector meson state or as 
decoupled partons if the px of the partons is high enough. The splitting 7 — > qq is called 
the anomalous component of the photon. The structure function of virtual photons has 
only recently been measured but by far not to the same precision as the proton structure 
function. However, it turns out that data are in good agreement with the parameterization 
of Schuler and Sjostrand (SaS) fL2| . The SaS parameterization offers a choice of Ql values 
at which the anomalous part becomes effective. We have studied these choices resulting 
in different magnitudes of the parton densities, and consequently of the cross sections. 
For the SaS parameterization we have used Ql as given by eq.(12) of ref. [|l^] (IP2 = 2). 
This choice is also suitable for a description of other hadronic final state properties (not 
considered in this paper), like energy flow, forward particle spectra and jet cross sections. 

In a previous paper n we have shown that the ansatz of Drees and Godbole [p3 



gave very similar results to those obtained from the SaS parameterization why we restrict 
ourselves to the SaS parton distributions in this study. 

The hadronic contribution to the virtual photon structure function decreases rapidly 
with increasing Q 2 , which means that the main contribution at large Q 2 comes from the 
anomalous piece in the photon splitting. This is completely calculable in pQCD and leads 



to an expected agreement between the parameterization of Gliick - Reya - Stratman JL1 
and that of Schuler - Sjostrand [|T!J but also with the simple ansatz of Drees - Godbole 
pi . However differences exist in the way the hadronic part of the structure function is 



matched to the pointlike part, which just reflects the theoretical uncertainty. 



2.2 Choice of Scale 

In leading order a s processes the renormalization scale fiR and factorization scale fip 
are not well defined which allows a number of reasonable choices. There are essentially 
two competing effects: a large scale suppresses a s (/z 2 ) but gives, on the other hand, an 
increased parton density, x/(x,/i 2 ), for a fixed small x value. The net effect depends on 
the details of the interaction and on the parton density parameterization. 

In previous papers (16], [L| we have tried different scales like /i 2 = and /i 2 = Q 2 +p?, 



and found that these choices gave similar results. 

However, in resolved virtual photon processes the choice of the scale fi 2 , at which the 
photon is probed, is severely restricted fl~7|l . In a partonic process a + b — > c + d, where 



a, b, c, d denote four-vectors and where parton a has the virtuality Q 2 , the transverse 
momentum p\ of parton c is given in the small angle limit (— t <C s) by: p\ = s(— t) /{s + 
Q 2 ), with s and t being the usual Mandelstam variables. In a t channel process the 
virtuality is given by fi 2 = —t. Thus we have:[] 



H = —t = Pt + Q ■ — < Q + p T 



s 



From eq.(0) we see, that the scale /i 2 is always larger than the transverse momentum 
squared of the hard partons and less than Q 2 +p^. In the following we will use fi 2 = Q 2 +p\ 
as scale for both resolved virtual photon processes and for direct photon processes. This 
choice of scale provides a smooth transition from the kinematic region of normal DIS into 
the range where resolved photons start contributing and to the photo-production region. 
The same scale has also been used in NLO calculations including resolved photons in deep 
inelastic scattering 



A basic test that the scale is reasonable is that the parton shower evolution scheme 
should not be able to produce partons with transverse momenta larger than those pro- 
duced by the matrix element for the hard scattering process. 



3 Forward Jets 

HERA has extended the available x region down to values below 10 -4 where new parton 
dynamics might show up. Based on calculations in the LLA of the BFKL kernel, the cross 
section for DIS events at low x and large Q 2 with a high p\ jet in the proton direction (a 
forward jet) || |3[ is expected to rise more rapidly with decreasing x than expected from 
DGLAP based calculations. New results from the HI B and ZEUS H experiments have 
recently been presented. The data can be described neither by conventional DIR Monte 
Carlo models nor by a NLO calculation, while comparisons to analytic calculations of the 
LLA BFKL mechanism has proven reasonable agreement. 

It should be kept in mind that both the NLO calculations and the BFKL based 
calculations are performed on the parton level whereas the data are on the level of hadrons. 

*We are grateful to T. Sjostrand for pointing out this simple explanation 
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Figure 2: The forward jet cross section as a function of x for pr jet > 3.5 GeV (a.) 
and Pt jet > 5 GeV (&.). Also shown are the RAPGAP predictions for the sum of 
direct and resolved processes (solid line) as well as the resolved photon contribution 
alone (dashed line). 

In Fig. H] the forward jet cross section as measured by the HI collaboration || is 
compared to the prediction of the RAPGAP Monte Carlo generator for both, resolved 
photon process alone (labeled RES) and for the sum of direct and resolved processes 
(labeled DIR+RES). The calculation is performed with a scale fi 2 = Q 2 + p\ used in the 
parton densities and a s for both direct and resolved photon processes. The measurement 
is well described by the Monte Carlo, including direct and resolved photon processes. 
The forward jet cross section as measured by the ZEUS experiment || can be equally well 
described by this Monte Carlo program using the same structure functions and parameter 



setting. 




Figure 3: The ratio R = q\jp\ °f the transverse momenta q T of partons from the 
initial state cascade to the transverse momentum p\ of the partons from the hard 
scattering process, the solid line corresponds to the scale fi 2 = Q 2 + p T , the dotted 
line to /i 2 = 4 ■ P% and the dashed line to fi 2 = p\. The distribution is normalized 
to the total number of entries. Please note the logarithmic scale on the y - axis. 

For the forward jet analysis it is important to check that the reconstructed jets really 
stem from the hard scattering and not from the initial state parton cascade. In Fig. [3] 
the ratio of the transverse momentum of any initial state parton q T to the transverse 
momentum p\ of the hard scattering process is shown in the 7*p CMS for events which 
satisfy the forward jet analysis criteria. The solid line corresponds to /i 2 = Q 2 + p^, the 
dotted line to fj, 2 = 4 -p\ and the dashed line to /i 2 = p\. One can see, that essentially all 
partons coming from the initial state cascade have transverse momenta smaller than the 
partons of the hard scattering px, which is expected in a DGLAP type evolution, where 
the transverse momenta are ordered in qx towards the hard scattering process. Thus 
we conclude that the scale /i 2 = Q 2 + p\ which has been used for this study fullfils the 



requirements of a DGLAP type initial state cascade. 
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Figure 4: Different contributions to the total cross section of resolved virtual photons 
within the cuts of the forward jet analysis H. In a. is shown the ratio R = — — — , 

' — ' <J res. tot 

i.e. the pointlike part (solid line) and the hadronic part (dashed line), respectively, 
of the resolved virtual photon cross section divided by the the total resolved photon 
cross section as a function of x. In b. is shown the ratio R = — — — as a function 

&res. tot 

of x for different subprocesses i: qg — > qg (dotted line), qq — ► qq and qq — > qq (solid 
line), gg — > (dashed line), qq — ► (dashed-dotted line) and gg — ► (solid line). 



In Fig. ^ the different contributions of the total resolved photon cross section are shown 
separately within the cuts of the forward jet analysis. From Fig. it is observed that 
the hadronic part of the virtual photon structure function as expected gives a negligible 



contribution to the measured cross section, since it dies off rapidly with increasing Q 2 . 
Fig. fp? shows that the subprocess q^g v — > qg contributes the most to the resolved photon 
cross section in the forward jet region (~ 60%) and that the subprocesses qq — > qq, 
qq — ► qq and gg — ► gg each give a contribution of the order of 10%. 

A small fraction of the DIS events, fulfilling the selection criteria for forward jets, 
actually contains two identified jets. Analytic calculations (in LLA) |HJ have been per- 
formed in the same kinematic region and with the same jet selection as defined for the 
one-jet sample. The predicted ratio varies from 3% to 6% as x increases from 0.5 • 10~ 3 
to x = 3 • 10~ 3 . Our previously reported prediction from the RAPGAP generator |IJ 
including both direct and resolved photon processes was that about 1 % of the total 
forward jet sample contains two forward jets. This is about a factor of 3 lower than 
the prediction from the BFKL calculations but a large part of this discrepancy could be 
due to hadronization effects which would reduce the prediction of the parton level BFKL 
calculation. 

Recently this ratio has been measured by the HI experiment || to be 1 %, in excellent 
agreement with the prediction of RAPGAP. This gives further confidence in the basic 
concept of resolved photons even at large Q 2 . 

The ZEUS collaboration has presented a measurement of the forward jet cross section 
as a function of p 2 T jQ 2 |2(J in the kinematic region Q 2 > 10 GeV 2 , y > 0.1, E e > > 10 
GeV, ri jet < 2.6, x jet > 0.036, E T jet > 5 GeV, p zbreit > GeV, 2.5 • 10~ 4 < x < 8 • 10~ 2 
but without implementing the DGLAP suppression cut 0.5 < E\jQ 2 < 2. The results 
are compared to the predictions from different Monte Carlo programs, and the conclusion 
is that only the RAPGAP DIS generator including interactions through resolved virtual 
photons can describe the data over the full range in E\jQ 2 . In Fig. |5|a the ZEUS results 
are shown together with the prediction of RAPGAP and in Fig. |5|& the contributions 
coming from the hadronic and the pointlike part of the virtual photon structure function 
are presented separately. 



4 Summary and Discussion 

Recent experimental data on forward jet production show deviations from traditional LO 
Monte Carlo models assuming directly interacting point-like photons. It is tempting to 
assume that the observed effects could be explained by BFKL dynamics. 

In the present study we have shown that the addition of resolved photon processes to 
the direct interactions in DIS leads to good agreement with the data. This agreement 
does not depend on any specific choice of scale or tuning of any other parameters in the 
RAPGAP generator. The best evidence of the universality of this approach is that, with 
the same parameter setting, it is possible to describe a wide range of other data like 
the transverse energy flow [21], transverse momentum spectra of single particles ||22|| , the 



(2+1) jet rate |23| and single inclusive jet cross sections ||, as we have shown in JT|. 

We have observed that the dominant contributions to the resolved photon processes 
come from order a 2 diagrams with the hard subprocess q^g p — > qg (see Fig. [[J). Since the 
partons which form the photon remnant per definition have smaller p? than the partons 
involved in the hard scattering, a situation with non q? ordering is created. 
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Figure 5: The cross section of forward jets as a function of E\jQ 2 . In a. is shown 
the measurement of ZEUS [20| and the prediction of RAPGAP. The solid line shows 
the sum of direct and resolved virtual photon contributions, whereas the dashed line 
shows the resolved photon contribution alone. In h. is shown the part of the cross 
section coming from the anomalous component (solid line) and the one coming from 
the hadronic component of the virtual photon separately. 



In the LO DIR model the ladder of gluon emissions is governed by DGLAP dynamics 
giving a strong ordering of q t for emissions between the photon and the proton vertex. The 
models describing resolved photon processes and BFKL dynamics are similar in the sense 
that both lead to a breaking of this ordering in q t . The BFKL picture, however, allows 
for complete dis-ordering in q t , while in the resolved photon case the DGLAP ladder is 
split into two shorter ladders, one from the hard subsystem to the proton vertex, and one 



to the photon vertex, each of them ordered in q t (see Fig. [T]). Only if the ladders are long 
enough to produce additional hard radiation it might be possible to separate resolved 
photon processes from processes governed by BFKL dynamics. Thus the resolved photon 
approach may be a "sufficiently good" approximation to an exact BFKL calculation and 
the two approaches may prove indistinguishable within the range of x accessible at HERA. 

It should be emphasized again that a NLO calculation assuming point-like virtual 
photons contains a significant part of what is attributed to the resolved structure of the 



virtual photon in the RES model [IS 
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